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Abstract 



These lectures aim at providing an introduction to the properties of gravitational 
waves and in particular to those gravitational waves that are expected as a consequence 
of perturbations of black holes and neutron stars. Imprinted in the gravitational ra- 
diation emitted by these objects is, in fact, a wealth of physical information. In the 
case of black holes, a detailed knowledge of the gravitational radiation emitted as a 
response to perturbations will reveal us important details about their mass and spin, 
but also about the fundamental properties of the event horizon. In the case of neutron 
stars, on the other hand, this information can provide a detailed map of their internal 
structure and tell us about the equation of state of matter at very high density, thus 
filling-in a gap in energies and densities that cannot be investigated by experiments in 
terrestrial laboratories. 
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1 A Brief Introduction to Gravitational Waves 

Being prepared for a Summer School, these lecture notes are meant principally as a first 
introduction to the rather vast and rapidly growing area of research that deals with gravita- 
tional waves. As a result, these notes will simply provide the basic concepts of this field of 
research and briefly review the most important results. The bibliography is not meant to be 
complete and the sparse references are used mostly as pointers to review articles where the 
topics of these lectures find a more detailed discussion and where a complete presentation 
of the literature can be found. In what follows I will assume that the reader is proficient 
in basic tensor algebra, is familiar with the fundamental concepts in General Relativity and 
has already encountered a black hole solution to Einstein equations. 

The lecture notes are organized in three main parts. In the first one, I rapidly introduce 
the basic properties of gravitational waves by deriving a wave solution of the Einstein equa- 
tions and discussing the sense of the transverse-traceless gauge. A brief discussion about the 
generation of gravitational waves will also be made. In the second part, I discuss in some 
detail the perturbation properties of a Schwarzschild black hole and the gravitational radi- 
ation that is expected as a result. Finally, the third part is devoted to presenting the basic 
elements of the theory of stellar perturbations and how stellar perturbations can be used to 
generate gravitational waves. The discussion will be restricted to spherical oscillations in the 
case of relativistic stars and to nonradial oscillations in the case of Newtonian stars. A brief 
classification of the modes of oscillations will also be presented together with an introduction 
to the onset of non-axisymmetric instabilities. 

These lectures notes are taken in part from the introductory course on General Relativity 
given to the graduate students at SISSA and is inspired by the corresponding discussions 
in a number of textbooks and in particular by those in ^ [2 El- A space-like signature 
(—,+,+, +) will be used, with Greek indices taken to run from to 3 and Latin indices 
from 1 to 3. Covariant derivatives are denoted with a semi-colon and partial derivatives 
with a comma. Tensors are indicated with bold symbols (i.e. A) and three-vectors with the 
standard arrow (i.e. A). Finally, I will adopt geometrized units in which G — c — 1. All of 
the figures presented here are original and have not been published elsewhere. 

1.1 Linearized Einstein Equations 

The starting point in discussing gravitational waves cannot but come from the Einstein field 
equations, expressing the close equivalence between matter-energy and curvature 



In the 10 linearly independent equations (Q), R^ v and R are the Ricci tensor and scalar, 
respectively, g^ u and G^ u are the metric and Einstein tensors, respectively, while T^ u is the 
stress-energy tensor of the matter in the spacetime considered. 

Looking at the Einstein equations (0) as a set of second-order partial differential equa- 
tions it is not easy to predict that there exist solutions behaving as waves. Indeed, and as it 
will become more apparent in this Section, the concept of gravitational waves as solutions of 
Einstein equations is valid only under some rather idealized assumptions such as: a vacuum 
and asymptotically flat spacetime and a linearized regime for the gravitational fields. If these 
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assumptions are removed, the definition of gravitational waves becomes much more difficult. 
In these cases, in fact, the full nonlinearity of the Einstein equations complicates the treat- 
ment considerably and solutions can be found only numerically. It should be noted, however, 
that in this respect gravitational waves are not peculiar. Any wave-like phenomenon, in fact, 
can be described in terms of exact "wave equations" only under very simplified assumptions 
such as those requiring an uniform "background" for the fields propagating as waves. 

These considerations suggest that the search for wave-like solutions to Einstein equations 
should be made in a spacetime with very modest curvature and with a metric line element 
which is that of flat spacetime but for small deviations on nonzero curvature, i.e. 

= V^ + V + ([VH > ( 2 ) 

where 

r},^ = diag(-l, 1, 1, 1) , (3) 
and the linearized regime is guaranteed by the fact that 

I VI < 1 ■ (4) 

Fortunately, the conditions expressed by equations (j2J) and (jU) are, at least in our Solar 
system, rather easy to reproduce and, in fact, the deviation away from flat spacetime that 
could be measured, for instance, on the surface of the Sun are 

I VI ~ I VI - ~ 10-6 • ( 5 ) 

Before writing the linearized version of the Einstein equations it is necessary to derive 
the linearized expression for the Christoffel symbols. In a coordinate basis (as the one will 
will assume hereafter), the general expression for the affine connection is 



r % = 7:9^(9^,(3 + 9/3u,a ~ 9af3,v) , (6) 



where the partial derivatives are readily calculated as 

+ h 

so that the linearized Christoffel symbols become 

= ^% + V,«-V M )- (8) 

Note that the operation of lowering and raising the indices in expression (jHJ) is not made 
through the metric tensors g^ u and g^ u but, rather, through the spacetime metric tensors 
rj^ and n >Ml/ . This is just the consequence of linearized approximation and, despite this, the 
spacetime is really curved! 
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Once the linearized Christoffel symbols have been computed, it is possible to derive the 
linearized expression for the Ricci tensor which takes the form 

p pa -pa 

= 2 V ' ua ~^^ l ' lJ 'A 4 " _ ^^v,a ~ h,fJ.v) ) (9) 

where 

h = h\ = rTh^ , (10) 
is the trace of the metric perturbations. The resulting Ricci scalar is then given by 

R^g^R^-^R^ . (11) 

Making now use of and (jllj) it is possible to rewrite the Einstein equations (0) in a 
linearized form as 

h, a % + K a ^ - h^ a - h, u - r]^(h aP f - h a a ) = 16vrT^ . (12) 

Although linearized, the Einstein equations ()12|) do not seem yet to suggest a wave-like 
behaviour. A good step in the direction of unveiling this behaviour can be made if a more 
compact notation is introduced and which makes use of "trace-free" tensors defined as 

h^v = h^ v - -rj^h , (13) 

where the "bar-operator" can be applied to any symmetric tensor so that, for instance, 
R^iv = G[iv and h^ u = h^ 1 . Using this notation, the linearized Einstein equations (fT2|) take 
the more compact form 

- W ~ V^apf + ha, \ = l&lT^ • (14) 

It is now straightforward to recognize that the first term on the right-hand-side of equation 
(J 14)1 is simply the Dalambertian (or wave) operator 

W = n V = +% + % + 9 2 z ) V , (15) 

where the last equality is valid for a Cartesian (t, x,y, z) coordinate system only. At this stage 
the gauge freedom inherent to General Relativity can (and should) be exploited to recast 
equations (|15|) in a more convenient form. A good way of exploiting this gauge freedom is 
by choosing the metric perturbations h^ u so as to eliminate the terms in (J 14)) that spoil the 
wave-like structure. Most notably, the metric perturbations can be selected so that 

^% = . (16) 

Making use of the gauge (|16|). which is also known as "Lorentz" (or Hilbert) gauge, the 
linearized field equations take the form 

□V = -16ttT^ . (17) 



1 Note that the "bar" operator can in principle be applied also to the trace so that h 
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Despite they are treated in a linearized regime and with a proper choice of variables and 
gauges, Einstein equations (|T7|) do not yet represent wave-like equations if matter is present 
(i.e. if 7^ 0). A further and final step needs therefore to be taken and this amounts to 
consider a spacetime devoid of matter, in which the Einstein equations can finally be written 

as 

□V = > ( 18 ) 

indicating that, in the Lorentz gauge, the "gravitational field" propagates in spacetime as a 
wave perturbing flat spacetime. 

Having recast the Einstein field equations in a wave-like form has brought us just half-way 
towards analysing the properties of these objects. More will be needed in order to discuss 
the nature and features of gravitational waves and this is what is presented in the following 
Section. 



1.2 A Wave Solution to Einstein Equations 

The simplest solution to the linearized Einstein equations (JTSj) is that of a plane wave of the 
type 

h^u = U {A^ v exp(iK a x a )} , (19) 

where A is the "amplitude tensor" and ft is a null four-vector, i.e. K a n a = 0. In such 
a solution, the plane wave (fT9Jl travels in the spatial direction k = (k x , n y , k z )/k° with 
frequency uj = k° = (k^kj) 1 ^ 2 . 

Note that the amplitude tensor A in the wave solution (jl9j) has in principle 16 — 6 = 10 
independent components. On the other hand, a number of considerations indicate that there 
are only two dynamical degrees of freedom in General Relativity. This "excess" of indepen- 
dent components can be explained simply. Firstly, A and n cannot be arbitrary if they have 
to describe a plane wave; as a result, an orthogonality condition between the two quantities 
will constrain four of the ten components of A (see condition (a) below). Secondly, while 
a global Lorentz gauge has been chosen [cf. equation (fTfij)]. this does not fix completely the 
coordinate system of a linearized theory. A residual ambiguity, in fact, is preserved through 
arbitrary "gauge changes", i.e. through infinitesimal coordinate transformations that are not 
entirely constrained, even if a global gauge has been selected. To better appreciate this, con- 
sider an infinitesimal coordinate transformation in terms of a small but otherwise arbitrary 
displacement four-vector £ 

x a ' = x a + £ a . (20) 

Applying this transformation to the linearized metric (J2J) generates a "new" metric tensor 
that at the lowest order is 

so that the "new" and "old" perturbations are related by the following expression 

-NEW , OLD , , 

hplvl — fytf, _ £(U,K ~ ■ (22) 

or, alternatively, by 

-NEW T OLD t /-a , , 

Vw 1 = iiv ~ S/i,^ ~ <iv,\l + ,a ■ V^J 
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Requiring now that the new coordinates satisfy the condition (|16p of the Lorentz gauge, 
forces the displacement vector to be solution of the sourceless wave equation 

C'% = • (24) 

As a result, the plane-wave vector with components 

i a = -iC a exp(iKpx p ) , (25) 

generates, through the four arbitrary constants C a , a gauge transformation that changes 
arbitrarily four more components of A. Effectively, therefore, A^ v has only 10 — 4 — 4 = 2 
linearly independent components, corresponding to the number of degrees of freedom in 
General Relativity PQ. 

Note that all this is very similar to what happens in classical electrodynamics, where 
the Maxwell equations are invariant under transformations of the vector potentials of the 
type — > v4y = A^ + \l/ >At , so that the corresponding electromagnetic tensor F , u , = 

Apiy — Ayi ^i = F°^ . Similarly, in a linearized theory of General Relativity, the gauge 

transformation (jUj) will preserve the components of the Riemann tensor, i.e. F^L, = ^a/w- 
in practice, it is often convenient to constrain the components of the amplitude tensor 
through the following conditions: 

(a) : Orthogonality Condition: four components of the amplitude tensor can be specified 

if A and k, are chosen to be orthogonal 

A^k v = . (26) 

(b) : Global Lorentz Frame: just like in Special Relativity, a global Lorentz frame relative 

to an observer with four-velocity u can be defined. In this case, three 2 components of 
the amplitude tensor can be specified after selecting a four-velocity u orthogonal to A 

A^ v u v = . (27) 

(c) : Infinitesimal Gauge Transformation: one final independent component in the 

amplitude tensor can be eliminated after selecting the infinitesimal displacement vector 
£ M = iC 1 exp(iK a x a ) so that 



A% = . (28) 

Consider now the constraint conditions (|26j) - (|28|) listed above as implemented in a ref- 
erence frame which is globally at rest, i.e. u a = (1,0,0,0). In this frame, the components 
of the wave vector do not appear directly, and the above conditions for the amplitude 
tensor can be written as 

2 Note that the conditions l|26|l fix three and not four components because one further constraint needs 
to be satisfied, i.e. n^A^u" = 0. 
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(a) : 

A^X = h ijtj = , (29) 

i.e. the spatial components of are divergence-free. 

(b) : 

A^u v = <^ hpo = , (30) 
i.e. only the spatial components of h^ u are nonzero. 

(c) : 

A^ = ^ h = h 3 j = 0, (31) 

i.e. the spatial components of h^ v are trace-free. Because of this, and only in this 
gauge, = 

Conditions (a), (b) and (c) define the so called "Transverse and Traceless" (TT) gauge and 
represent the standard gauge for the analysis of gravitational waves. 

An obvious question that might emerge at this point is about the generality of the TT 
gauge. A simple answer to this question can be provided by reminding that any linear 
gravitational wave can, just like any electromagnetic wave, be decomposed in the linear 
superposition of planar waves. Because all of the conditions (jSUjl - piJl are linear in h^, any 
of the composing planar waves can be chosen to satisfy (j3TH) - (}3~Tj) . which, as a result, are 
satisfied also by the original gravitational wave. Indeed, all of what just stated is contained in 
a theorem establishing that: once a global Lorentz frame has been chosen in which u a = S a Q , 
it is then always possible to find a gauge in which the conditions fl3T H) - (j3Tj) are satisfied. 



1.3 Making Sense of the TT Gauge 

As introduced so far, the TT gauge might appear rather abstract and not particularly inter- 
esting. Quite the opposite, the TT gauge introduces a number of important advantages and 
simplifications in the study of gravitational waves. The most important of these is that, in 
this gauge, the only nonzero components of the Riemann tensor are 

RjOkO — RojOk — —RjOOk — —RojkO ■ (32) 

Since, however, 

RjOkO = -gfyfc.oo ' ( 33 ) 

the use of the TT gauge indicates that a travelling gravitational wave with periodic time 
behaviour h^ k oc exp(iut) can be associated to a local oscillation of the spacetime, i.e. 

^JJoo ~ -uj 2 exp(iut) ~ R j0k0 , and R j0k0 = -o; 2 /iJ^ . (34) 

To better appreciate the effects of the propagation of a gravitational wave, it is useful to 
consider the separation between two neighbouring particles A and Bona geodesic motion and 
how this separation changes in the presence of an incident gravitational wave (see Figure Q). 
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For this purpose, let us introduce a coordinate system x a in the neighbourhood of particle 
A so that along the worldline of the particle A the line element will have the form 

ds 2 = -dr 2 + ,) ;j d.r'd.r j + 0(\x j \ 2 )dx & dx? . (35) 

The arrival of a gravitational wave will perturb the geodesic motion of the two particles 
and produce a nonzero contribution to the geodesic-deviation equation. I remind that the 
changes in the separation four-vector V between two geodesic trajectories with tangent 
four-vector u are expressed through the geodesic-deviation equation 

v?u^V a . M = -R^y^u 5 , (36) 

or, equivalently, as 

7 p(D 2 V a \ _ p(d 2 V a dV a dVn p s 

{-^- ) = ^ + r *-*r-*r) = ~ R ^ u V u ■ (37) 




Figure 1: Schematic diagram for the changes in the separation vector between two particles A and B 
moving along geodesic trajectories produced by the interaction with a gravitational wave propagating along 
the direction k. 



Indicating now with = x J B — = x J B the components of the separation three-vector 
in the positions of the two particles, the geodesic-deviation equation (|37j) can be written as 

D 2 x~ j 

= -*Wb • (38) 

A first simplification to these equations comes from the fact that around the particle A the 
affine connections vanish (i.e. F{~ = 0) and the covariant derivative in ()38|) can be replaced 
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by an ordinary total derivative. Furthermore, because in the TT gauge the coordinate system 
x a moves together with the particle A, the proper and the coordinate time coincide at first 



order in the metric perturbation [i.e. r 
becomes 

d 2 xl 



t at 0(h )]. As a result, equation (JHS|) effectively 



dt 2 



'd 2 h 



dt 2 



X: 



(39) 



and has solution 



xUt)=xU0) 



1 , TT , . 

¥ + 2% (*) 



(40) 



Equation (|4U|) has a straightforward interpretation and indicates that, in the reference frame 

TT 

comoving with A, the particle B is seen oscillating with an amplitude proportional to h-~ . 

Note that because these are transverse waves, they will produce a local deformation of 
the spacetime only in the plane orthogonal to their direction of propagation. As a result, 

k), then /^V B (0) oc 



if the two particles lay along the direction of propagation (i.e. if n 



^4(0) 



and no oscillation will be recorded by A [cf. equation (f2T)|) ] 






t=P/4 



t=P/2 



t=3P/4 



t=P 







Figure 2: Schematic deformations produced on a ring of freely-falling particles by gravitational waves that 
are linear polarized in the "+" ("plus") and "x" ("cross") modes. The continuous lines and the dark filled 
dots show the positions of the particles at different times, while the dashed lines and the open dots show the 
unperturbed positions. 



Let us now consider a concrete example and in particular a planar gravitational wave 
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propagating in the positive ^-direction. In this case 

TT TT 

K x = -h yy =&{A+exp[-iu(t-z)]} , (41) 

TT TT 

K y = h yx =M{A x ex V [-iu{t-z))} , (42) 

where A + and A x represent the two independent modes of polarization. 
As in classical electromagnetism, in fact, it is possible to decompose a gravitational wave in 
two linearly polarized plane waves or in two circularly polarized ones. In the first case, and 
for a gravitational wave propagating in the z-direction, the polarization tensors + ("plus") 
and x ( "cross" ) are defined as 

e + = e x ®e x -e y ®e y , (43) 

e x = e x ® e x + e y ® e y . (44) 

The deformations that are associated with these two modes of linear polarization are 
shown in Figure El where the positions of a ring of freely-falling particles are schematically 
represented at different fractions of an oscillation period. Note that the two linear polariza- 
tion modes are simply rotated of 7r/4. 

In a similar way, it is possible to define two tensors describing the two states of circu- 
lar polarization and indicate with e R the circular polarization that rotates clockwise (see 
Figure OJ) 

and with e L the circular polarization that rotates counter-clockwise (see Figure |3J) 

e L = . (46) 

The deformations that are associated to these two modes of circular polarization are shown 
in Figure El 

1.4 Generation of Gravitational Waves 

In what follows I briefly discuss the amounts of energy carried by gravitational waves and 
provide simple expressions to estimate the gravitational radiation luminosity of potential 
sources. Despite the estimates made here come from analogies with electromagnetism and 
are based on a Newtonian description of gravity, they provide a reasonable approximation 
to more accurate expressions from which they differ for factors of O (few). Note also that if 
obtaining such a level of accuracy requires only a small effort, reaching the accuracy necessary 
for a realistic detection of gravitational waves is far more difficult and often imposes the use 
of numerical relativity calculations on modern supercomputers. 

In classical electrodynamics, the energy emitted per unit time by an oscillating electric 
dipole d is easily estimated to be 

_ (energy emitted) 2 22 2 •■ 2 

-^electric dip. = 7 r— r — ~q a — -{d) , (4 ( ) 

unit time 3 3 
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t=P/4 



t=P/2 



t=3P/4 



t=P 





Figure 3: Schematic deformations produced on a ring of freely- falling particles by gravitational waves that 
are circularly polarized in the R (clockwise) and L (counter-clockwise) modes. The continuous lines and the 
dark filled dots show the positions of the particles at different times, while the dashed lines and the open 
dots show the unperturbed positions. 



where 

d = qx , and d = qx , (48) 

with q being the electrical charge and the number of "dots" counting the order of the total 
time derivative. Equally simple is to calculate the corresponding luminosity in gravitational 
waves produced by an oscillating mass-dipole. In the case of a system of N point-like particles 
of mass rriA (A = 1, 2, . . . , N), in fact, the total mass-dipole and its first time derivative are 



N 



d = ^2 m A ^ A ' 

A=l 



(49) 



and 



N 



d = ^2 m AXA = P , (50) 

A=l 

respectively. However, the requirement that the system conserves its total linear momentum 



d = p = 0, (51) 

forces to conclude that L mass di po ie — 0, i.e. that there is no mass-dipole radiation in General 
Relativity (This is equivalent to the impossibility of having electromagnetic radiation from an 
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oscillating electric monopole.). Next, consider the electromagnetic energy emission produced 
by an oscillating electric quadrupole. In classical electrodynamics, this energy loss is given 
by 

-^electric quad. = — (Q) 2 = l^(QjkQjk) 2 J (52) 

where 



N 



Qjk = 1a 



A=l 



{x A ) j (x A ) k - ^S jk (x A )i( y x A y 



(53) 



is the electric quadrupole for a distribution of N charges (qi, q 2 , . . . , (Mr). 

In close analogy with expression (jo^j) . the energy loss per unit time due to an oscillating 

mass quadrupole is calculated to be 



] Q ... ] Q 

Lmass quadrupole = ~Z ~ g ( I ) T — 5 \ ^~ j'fe jk) i 

DC DC 



where fj k is the trace-less mass quadrupole (or "reduced" mass quadrupole), defined as 

N 



i 

{x A )j{x A ) k - -S jk (x A )i(x A y 



A=l 

p (xjX k - ^5 jk XiX^j dV , (55) 

and the brackets ( ) indicate a time average [Clearly, the second expression in ()55j) refers to 
a continuous distribution of particles with rest-mass density p.}. 

A crude estimate of the third derivative of the mass quadrupole of the system is given by 

£ (mass of the system in motion) x (size of the system) 2 MR 2 Mv 2 
3 (time scale) 3 r 3 r ' 

so that 

GMv 2 



7 mass quadrupole g ' 



(57) 

C J T 

Although extremely simplified, expressions (J5T)j) and (jBTf) contain the two most important 
pieces of information about the generation of gravitational waves. The first one is that the 
energy emission in gravitational waves is severely suppressed by the coefficient G/c 5 ~ 10~ 59 
or, stated differently, that the conversion of any type of energy into gravitational waves is, 
in general, not efficient. The second one is about the time variation of the mass quadrupole, 
which can become considerable only for very large masses moving at relativistic speeds. 
Clearly, these conditions cannot be reached by sources in terrestrial laboratories, but can 
be easily met by astrophysical compact objects, which therefore become the most promising 
sources of gravitational radiation. 



1.4.1 Astrophysical Sources of Gravitational Waves 



The research area that is involved with modelling the astrophysical sources of gravitational 
waves is quite vast and is increasing steadily as gravitational detectors become more and 
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more sensitive, and as detectors of new generation are becoming operative. Having a detailed 
discussion of the multiple aspects of this line of research, which encompasses numerical and 
perturbative techniques, is beyond the scope of these lectures. However, the interested reader 
will find more detailed discussions in the following review articles [HE]. 
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2 Gravitational Waves From Perturbed Black Holes 

This Section is dedicated to the analysis of the perturbations that characterise a black hole 
and more precisely a nonrotating (or Schwarzschild) black hole. Hereafter I will assume that 
the reader is familiar with the basic properties of such a black hole as a static solution to 
Einstein equations in a spherically symmetric and vacuum spacetime. More precisely, I will 
consider as known the fundamental concepts about physical and coordinate singularities, 
about the existence of an event horizon as well as Birkhoff's theorem on the uniqueness of 
Schwarzschild solution. All of these concepts are clearly discussed in [TJ |2J E] . 

Before discussing in detail black hole perturbations, one might wonder why black hole 
perturbations are interesting at all. Indeed, there are a number of good reasons why it 
is interesting and important to consider black hole perturbations. Firstly, the presence of 
perturbations can break the static properties of a black hole spacetime and be therefore at 
the origin of gravitational wave emission. Secondly, the gravitational waves emitted by a 
black hole carry information about its properties such as mass, spin and charge. Finally, by 
investigating the response of black holes to perturbations it is possible to deduce important 
conclusions on the stability of these objects (being a solution of Einstein equations, in fact, 
is just a sufficient condition for stability). 

Because this is such an important area of research, some of the main results date back 
to the first studies made by Regge and Wheeler and the subsequent developments that 
took place in the 70's El- 



2.1 Linear Perturbations of Black Holes 

The starting point in the analysis of black hole perturbations is, of course, the unperturbed 
solution represented by a Schwarzschild black hole with line element 



ds 2 = g dx»dx v = - ( 1 - — ) dt 2 + ( 1 - — ) dr 2 + r W , (58) 



2M\ , „ / 2M\~ 1 



and where g v represents the metric tensor of the static background spacetime. Because the 
latter is assumed to be vacuum, the Einstein equations assume a more compact form and 
can be written as 

^ = 0, (59) 



where R^ u is the Ricci tensor built with the background metric g . If small perturbations 
h^ u are now introduced, the resulting metric will be 

9/w = V + K» > ( 60 ) 
where, again, the perturbations considered are much smaller than the background, i.e. 

IVl/lvl « 1 ■ (61) 

Just as for the static background, the behaviour of the perturbed spacetime will be 
expressed by the Einstein equations that, using the same notation as in ((HSJ), can be written 

as 

Rp, = , (62) 
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where R^ u = R^ u {g^ v )- At first order in the perturbations, the linearity can be exploited to 
break up the Einstein equations (|62|) as 

Rp, + 8R P v = , (63) 

where SR^ U = R^, v {h^ v ). Using now equations the field equations reduce to 

SRp, = . (64) 

A different way of writing equations (jSljl is in terms of the Christoffel symbols and more 
precisely as 

Mr = -K^ + = > ( 65 ) 

where the perturbed Christoffel symbols are defined as 

^{L = \g aP ( h fia,v + Ka,n ~ V,a) • ( 66 ) 

Before proceeding further I should comment on some of the properties of the metric per- 
turbations hfj, v . An important constraint is posed by Birkhoff 's theorem, which states that 
the Schwarzschild solution is the only spherically symmetric, asymptotically flat solution 
of Einstein equations in vacuum even if the spacetime is not static P3 121 El- As a result, 
nonrotating black holes can only be perturbed by nonradial perturbations and this forces to 
consider perturbations with complete angular dependence, i.e. = (t,r,9,(p). Han- 
dling a generic angular dependence can be complicated, but the mathematical treatment 
can be simplified if the tensor perturbations h^ u are written in a separable form i.e. as the 
product of four parts each being a function of one coordinate only. 

In the case of a scalar function depending on the spatial coordinates only, it is well-known 
that this can be done after expanding it in a series of spherical harmonic functions 

f(r,9,4 > ) = J2^m(r)Y lm (9,4 > ) . (67) 

l,m 

In a similar way, in the case of a vector, the separability is achieved through an expansion 
in a series of vector spherical harmonics [TU] 

V-(r,M) = ^a /m (r)[^(^,0)]« + ^^ m (r)[F^(^,0)] a , (68) 

l,m l,m 

where Y^(9, <fi) and Y^(9, 0) are vector spherical harmonics of magnetic (B) and electric 
(E) type, respectively. It will not therefore surprise that a series expansion of the type (Jo"7j) 
and (J68)) can be used also for a rank-2 symmetric tensor which can then be expanded in a 
series of tensor spherical harmonics 



T^(t, r, 9, 4>) = aim(t, r) [A^(9, 4>)], u + ^ b lm (t, r) [b£{0, (/>)] , (69) 

l,m l,m 
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where attention has been paid to the fact that, in general, a rank-2 symmetric tensor can 
be expanded in terms of tensor spherical harmonics that behave differently under parity 
transformation, i.e. (Af^)^ and (Bf^)^. 

More specifically, if V is the parity operator, that is an operator producing a parity 
transformation on a rank-2 symmetric tensor F^ v 

v([F lm (0,<f>))) — [^(^-0,^ + 0)1 , (70) 



J fllS 



the tensor spherical harmonics can then be classified according to their behaviour "under 
parity change". In practice, are referred to as odd or axial (or sometimes toroidal) those tensor 
harmonics for which V (F^ v ) = F^ u = (— 1)' +1 F^. Similarly, are referred to as even or polar 
(or sometimes spheroidal) those tensor harmonics for which V (F^ u ) = F^ v = (—I) 1 F^ u . 

This classification of the tensor spherical harmonics is reflected also on the metric per- 
turbations that, as a result, are classified as "odd" and "even-parity" respectively. However, 
before discussing the specific forms that the perturbations assume in these cases, it is con- 
venient to introduce the so-called "3 + 1" decomposition, in which the spacetime is "sliced" 
into a family of space-like spatial hypersurfaces parametrized by a t = const, coordinate. 
With this decomposition the line element takes the form 

ds 2 = -a 2 dt 2 + - fij (dx i + p i dt)(dx j + p j dt) 



= -(a 2 - p j pj)dt 2 + 2p i dx i dt + j ij dx i dx j , (71) 

where a is the "lapse" function (expressing the rate at which clocks on different hypersu- 
faces tick) and f3 l are the components of the "shift" vector (relating coordinate changes on 
two different hypersurfaces). With this decomposition, the metric perturbations are then 
expressed in terms of a purely time part (Z^oo); of a purely spatial part (hij), and of a mixed 
time-space part (h^) 

h - - ( t h C ) • < 72 > 

In what follows I discuss the basic expressions of ^ioo, h^j and hij for the two classes of 
perturbations as well as the equations they satisfy. 



2.2 Odd-parity Perturbations: the Regge- Wheeler Equation 

I first consider the parts of the metric perturbations that are of "odd-parity" [i.e. the first 
terms in the decomposition given in equation (J69j) ]. In this case, it is customary to introduce 
the unknown functions h (t,r), hi(t,r) and h 2 (t,r) so that the components of (f72|) can be 
written as 

h 00 = , (73) 



'0; 



h (t,r) 



, — 7—^ d^Yim, sin 9 V] d e Yi m 



Lm 



Lm 



h^ = h 1 (t,r)(e 1 ) ij + h 2 (t,r)(e 2 ) ij , 



(74) 
(75) 
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where (ei j2 )y = J2i m [(^i,2)ij]/ m an d I wm hereafter omit the l,m indices and the sum over 
them to maintain the expressions compact. The tensor spherical harmonics (ej 2)^ in 1)73)1- 
(1751) have rather lengthy but otherwise straightforward expressions which are given by 



/ 



ei. 







d^Yim sm9d 6 Yi m x 



sin 9 



sin # 



d<f>Yi m 



\ sin 6><9 e Y/ 



(76) 



and 



/ 



e 2 . 



\ 



sin 6 1 



(<9^ - cot 0fy) F Zm 



sin 2 



-di-cosOdo-smBdl 



Y, 



In 1 



\ 



1 

2 



-dl- cos edo-smOdl 



Y, 



1 in 



[sin 6><9^, - cos 98$] Y { 



lm 



• 2n~<P u o ~ L oca <fj - «« , 

Sill 6 v J 

(77) 

The Einstein equations with the metric perturbations (|73j ) -()77 p can be simplified if suit- 
able gauge conditions are chosen. I remind, in fact, that because of the linearized approach, 
any infinitesimal coordinate transformation will lead to new metric perturbations that are 
determined after the specification of the suitable conditions for the displacement four-vector 
£ M [cf. equations ()20j ) -()22j ) ]. While these conditions are totally arbitrary, it is convenient to 
choose those producing a simplification of the equations and, in the case of odd perturbations, 
the choice usually made is that 

h 2 (t,r) = 0. (78) 



In this gauge, which is usually referred to as the "Regge- Wheeler" gauge 
metric perturbations assume the simplified form 



/i ax 

"'flU 



( 









hi 










h \ 


/ 



sin 9dgP[ (cos 9) e 



imcj> 



the odd-parity 



(79) 



where Pi (cos 9) is the Legendre polynomial of order I. Besides being simpler, Einstein 
equations in this gauge are independent of m in the sense that the final result will not 
depend on the specific value chosen for m, which can therefore set to be zero. As a result, 
the Einstein equations for the perturbed metric ()79*j) lead to the following system of equations 



+ 1 



d 2 Q d 2 Q 
dt 2 dr 2 

dh d 

-dt- = dr-* inQ) 



2M 



1(1 + 1) 6M 



Q = 0, 



11) 



17 



where 



and 



r \ r 



r* = r 



2Mln 



2M 



- 1 



(82) 



13) 



is the "tortoise coordinate". Because r* — > r for r — > oo and r* — > — oo as r — > 2M + , the 
tortoise coordinate is particularly suited to study the propagation of perturbations near the 
black hole event horizon which, in this coordinate system, is placed at — oo and therefore 
does not suffer from coordinate singularities. 

A convenient way of looking at the "Regge- Wheeler equation" ()80|) is that of considering 
it as a wave equation in a scattering potential barrier V(r), where 



V(r) 



2M 



1(1 + 1) 



+ p 



with 



V 



6M 



(84) 



The potential (}84|) is also referred to as the "Regge- Wheeler potential" and has a maximum 
just outside the event horizon, at r ~ 3.3M in Schwarzschild coordinates, as shown in the 
schematic representation in Figure 0] 




Figure 4: Schematic representation of the Regge- Wheeler potential V [cf. eq. I)84|l ] and of its effect on a 
incident perturbation that is partially transmitted and partially reflected. 
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In this view, the Regge- Wheeler equation shares all of the well-known properties of a 
wave equation in a scattering potential and the numerous results that have been found for 
this type of equation (cf. Schrodinger equation) can also be applied to the propagation of 
perturbations in the spacetime of Schwarzschild black hole (see jTT] for a detailed discussion). 

As an example, a metric perturbations reaching the black hole from spatial infinity can 
be regarded as a wave packet that will scatter against the potential barrier V. As in quantum 
mechanics, not all of the wave packet will be transmitted through the potential and some 
of it, depending on the properties of the packet itself, will be reflected and reach again 
spatial infinity (see also Figure H}. This is different from what happens, for instance, with a 
spherical massive shell falling radially onto the black hole. These radical differences underline 
the importance of a perturbative analysis of black hole spacetimes. 

Another interesting aspect of the Regge- Wheeler equation is that it holds in a similar form 
also for scalar and vector perturbations, with the only difference appearing in the effective 
potential (|84|). where p = 2M/r 3 for scalar perturbations and p = for vector ones 



2.3 Even-parity Perturbations: the Zerilli Equation 

Next, I consider metric perturbation that are "even-parity" (or polar). The mathematical ap- 
proach is similar to the one followed for the odd-parity perturbations and also in this case it is 
useful to introduce a number of unknown functions h (t, r), h%(t, r), H (t, r), H^t, r), H 2 (t, r), 
K (t, r) and G(t, r) so that the perturbed metric functions can be written as 



hoi 



hi 



ij 



2M 



1/2 



H {t,r)Y lm , 



[H{Yi m , h doYi m , hod^Yim 
H 2 



hl{fl K, + l-2M/r VVij 



ij \ j % 3 



(85) 
(86) 

(87) 



where the tensor spherical harmonics (/ 1 _ 4 ) ij - have the forms 

/ d$Yi m dfkYim \ 



h 



d 6 Y lm 
V d^Y lm ) 



h 



'j 



/4 o o\ 

ooo 
\ o oo/ 



(89) 
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(90) 



and 



( \ 

d 2 e Y lm (d 2 0r/) - cot 9d,p)Y lm 

V {d 2 H - cot 6d^) Y lm (d 2 - sin 9 cos Odg) Y lm ) 



(91) 



Also for even-parity perturbations, the gauge freedom can be exploited and in particular 
a gauge can be chosen in which 



G = h = hi = . 



(92) 



As a result of this gauge choice, the polar metric perturbations assume the more compact 
form 



ol 



h p ' 



( H (l- 2M/r) Hi 

#i H 2 (1 - 2M/ry 1 

r 2 K 



\ 















,2-2, 



\ 



Pi (cos 6) e i 



(93) 



r 2 sin 2 9K j 



Writing now out the Einstein equations for the perturbed metric (|9*3^) leads to the following 
equation 

d 2 z d 2 z 



dt 2 Or 2 



VZ = 



(94) 



which is also known as the "Zerilli equation" . The explicit form of the "Zerilli function" is 
rather involved but can be expressed, independently of the gauge chosen, as 



Z = 



Are~ AX k 2 + 1(1 + l)rki 
1(1 + 1) -2 + 6M/r 



(95) 



where e A = 1 — 2M/r and the functions k±, k 2 , k 3 , fc 4 are introduced in place of G, h%, K, H 2 
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and are denned through the relations 



G = fc. 



(96) 



K 



Ho 



-2A 



2e~ 2X k 2 + r 



r 2 d h 
Or 

dki 
dr 



2k, 



1 +r— h-e- x — 
or J or 



2 -A ^3 



r e 



dr 



2e' x k 4 



(97) 
(98) 
(99) 



Note that as for odd-parity perturbations, the Einstein equations can be recast in the 
form of a wave equation in a scattering potential barrier V, defined as 



V= 1 



2M 



2q (q + 1) r 3 + 6q 2 Mr 2 + 18qM 2 r + 18M 3 
r 3 (qr + 3M) 2 



(100) 



where q = (l-l)(l + 2) /2. 

Interestingly, the Regge- Wheeler an Zerilli equations (JHUj) and are closely related 
and it is possible to transform the first one for axial modes into the second one for polar 
modes via suitable differential operators [TT] . 



2.4 QNMs of Black Holes 

Since equations (|50jl and ()94|) describe the response of the black hole to external pertur- 
bations, they are basically telling us about the vibrational modes of such a spacetime. In 
particular, if a harmonic time dependence is introduced for the perturbations in equations 
(|8U|l and i.e. if Q, Z ~ exp(ia; n t) where u;„ is the oscillation frequency of the n-th mode 
and is a complex number of the type 

u n = u r>n + iu iin , with n = 0,l,2,... , (101) 

it is then possible to define the Quasi-Normal Modes (QNMs) of the black hole as the 
solutions of equations 

d 2 t Q+[io 2 -V]Q = , (102) 



dl,Z + 



uo 2 -V 



, (103) 



that satisfy a pure outgoing-wave boundary condition at spatial infinity and a pure ingoing- 
wave boundary condition at the event horizon, i.e. 

Q, Z ~ exp(zu;r*) for r* — > — oo , and Q, Z ~ exp(— iur*) for r* — ► oo . (104) 
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2.4.1 A Summary of Main Results 

The literature on the solution of the Regge- Wheeler and of the Zerilli equations as well as 
on the determination of the perturbation spectrum of black holes is vast. The interested 
reader will find detailed discussions on these topics and the relevant references in the review 
by Kokkotas and Schmidt and in the one by Ferrari m what follows, however, I 
briefly summarize what could be considered the main results in the solution of the eigenvalue 
problem for the QNMs of a Schwarzschild black hole: 

• All the QNMs of Schwarzschild black hole have positive imaginary parts and represent 
therefore damped modes. As a result, a Schwarzschild black hole is linearly stable 
against perturbations. 

• The damping time of these perturbations depends linearly on the mass of the black 
hole (i.e. u n ~ 1/M) and is shorter for higher-order modes (i.e. u i;n+ i > to ijV ). As a 
result, the detection of gravitational waves emitted from a perturbed black hole could 
provide a direct measure of its mass. 

• The excitation of a black hole and the consequent emission of gravitational radiation is 
referred to as black hole "ringing" . The amplitudes of the gravitational waves emitted 
during the black hole ringing decay in time and the late-time behaviour (or tail of the 
ringing) is such that the amplitude evolution can be described in terms of a power law, 
whose envelope represents the superposition of the various QNMs. 

• The QNMs in black holes are isospectral, i.e. axial and polar perturbations have the 
same complex eigenfrequencies so that the real and imaginary parts of the spectrum 
are identical. This is simply due to the uniqueness in which a black hole can react to 
a perturbation. This is not true for relativistic stars. 

• The fundamental frequencies of oscillations u r ^ n /27r have been computed by a number 
of authors and are now known up to very large mode numbers. Reported in the table 
below are the frequencies of the first four modes, together with the values of their 
decaying times ra (i.e. I/uj^. The data refers to modes with I = 2, 3 and have been 
computed for M = M Q . 
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5.5344 x 10~ 2 
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5.3135 x 10- 2 
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11.203 


1.7983 x 10~ 2 
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18.833 


1.7510 x 10" 2 
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2 


9.7291 


1.0298 x 10~ 2 
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2 


17.834 


1.0281 x 10~ 2 


2 


3 


8.1264 


6.9856 x lO" 3 


3 


3 


16.551 


7.1354 x lO- 3 



• For any value of the harmonic index I, the real part of the frequency oo r>n (l) approaches 
a nonzero limiting value as the mode number n increases, while the imaginary part 
increases linearly as ~ n/4 (i.e. higher modes have shorter decaying timescales). 



Gravitational Waves from Perturbed Black Holes and Neutron Stars 

Most of what has been presented in this Section for a Schwarzschild black hole can be 
formulated also for a rotating (or Kerr) black hole. In this case, however, the mathe- 
matical apparatus is more involved (the potential is, for instance, complex) and some 
new features, such as the super-radiance (i.e. the amplified scattering of electromag- 
netic waves) can take place • More on this can be found in ^21 E] 
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3 Gravitational Waves from Perturbed Stars 

Black holes are often considered the most "extreme" objects predicted by General Relativity 
in the sense of being the objects with the most intense gravitational fields. This is of 
course correct, but it is worth reminding that relativistic stars (such as neutron stars or 
strange stars) are comparable to black holes in this respect and certainly the most relativistic 
"astrophysical" objects in the Universe. An important way in which they differ from black 
holes, however, is in possessing a surface and an interior structure which is not causally 
disconnected from the exterior spacetime. A rough estimate of how relativistic (and therefore 
"extreme") a compact star can be, is provided by its compactness, i.e. by the ratio between 
its gravitational mass M and its radius R. Shown in the table below is this quantity for a 
number of different objects 





M/R 


solar-type star 
white dwarf 
neutron star 
black hole 


~ lir 6 

~ 1(T 4 
0.1-0.3 
0.5 



Clearly, while an ordinary star cannot be considered a compact object, the table above 
shows that just a small difference in compactness distinguishes a neutron star from a black 
hole and this justifies the interest that relativists have reserved to this type of objects. In 
order to have a more concrete example of how "compact" a relativistic star can be, I have 
shown in Figure a schematic representation of the spatial dimensions of a neutron star 
superposing its surface to a local map of town for comparison. One should therefore imagine 
that an amount of matter a million times larger than the one contained in the Earth can be 
concentrated in an object with a radius of just ten kilometres. 

Just like black holes, compact stars can respond to perturbations and emit gravitational 
waves. Quite differently from black holes, however, the gravitational radiation coming from 
perturbed relativistic stars is extremely rich of details, many of which carry important phys- 
ical information. Imprinted in these gravitational waves, in fact, is a detailed map of the 
internal structure of the emitting stars, which can be used to deduce the properties of matter 
at conditions that cannot be investigated by experiments in terrestrial laboratories. 

Despite these exciting perspectives, the issue of the detectability of the gravitational 
radiation from perturbed relativistic stars is still basically unsettled. This is largely due to 
our ignorance about the precise physical conditions leading to a perturbed relativistic star. A 
simple example in this sense is offered by a protoneutron star formed after the gravitational 
collapse in a supernova explosion. While it is generally expected that the newly born neutron 
star will pulsate wildly during the first few seconds following the collapse, how much energy 
will be transferred to the pulsation and subsequently radiated through the oscillation modes 
is unknown. The only realistic way of overcoming this ignorance is to perform detailed, fully 
relativistic calculations and deduce from them how large the perturbations will be. While the 
recently developed numerical codes will soon be able to provide some quantitative answer in 
this respect JTHJ [T7| [TBI EH 1201 EE] , at present one can simply argue that a) the energy stored 
in the pulsation can potentially be of the same order as the kinetic energy of the system; b) 



24 



Gravitational Waves from Perturbed Black Holes and Neutron Stars 





Figure 5: Schematic representation of the spatial dimensions of a neutron star whose surface is shown as 
an open circle and is superposed over a local map of town for comparison. (This figure has been elaborated 
from one available on the AS-ICTP website.) 



the oscillations will be damped mainly through the emission of gravitational waves so that 
the released energy could be considerable. 

Under these assumptions, the effective gravitational wave amplitude h for a star oscillat- 
ing in its fundamental mode of oscillation can be estimated simply. For weak gravitational 
waves, in fact, the gravitational wave luminosity (i.e. the rate of energy loss to gravitational 
waves) at a distance d, can be written to be 



(,L 1:1 4nd 2 ( ) \h\ 2 , (105) 



dt t \16nG 

where E is the total energy lost over the time r. If the star is oscillating at a frequency /, 
then h ~ 2nfh and expression ()105|) can be rewritten as 



21 / E\ 1/2 ( 1 ms\ 1/2 ( 1 kHz\ / 50 kpc\ . . 

ft - 12x10 (l) (— ) {—){-?-)• (106) 

where E = 8.2 x 10~ 8 M Q c 2 is the energy lost in gravitational waves as estimated through 
recent relativistic calculations [XJJ]- Note that the probability of detecting a source can be 
increased if suitable data analysis techniques, such as "matched filtering", are used jUEj- 
In this case, it is possible to estimate the "effective" gravitational wave amplitude h e g to be 
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h cS ~ h^/Jr, so that ()106|) becomes 



^ I, x 10 - (^V'V^V'V^ ) . (107) 



,EJ \ f J \ d 

The distance scale d used in expressions (I106|) and (jl(J7|) is that to the supernova SN1987A 
and the number of events in the corresponding volume is one every 10-20 years. Clearly, this 
event rate is too small for being of interest and it is therefore necessary to consider a volume 
much larger, such as the one comprising the Virgo cluster, to reach an event rate of a few 
per year. In this case, it is possible to consider the problem of the detection from a different 
point of view and rather calculate what is the energy E necessary to obtain an effective wave 
amplitude h c s ~ 10~ 21 for a source at a distance d = 20 Mpc. Using expression ()107J) . the 
answer is E w 0.01 M . While these estimates may appear optimistic, they cannot be ruled 
out and the pay-offs of a potential detection would be so great to justify the intense research 
this field is experiencing. 

To better appreciate how the detection of gravitational waves from pulsating relativistic 
stars can be used to derive information about the physical properties of the source it is 
necessary to consider in more detail the perturbation theory of fluid stars, both in relativistic 
and in Newtonian gravity. This is the scope of the following Sections, where I will first 
discuss the spherical oscillations of relativistic stars and then the nonradial oscillations of 
Newtonian stars. A brief classification of the oscillation modes and an introduction to the 
onset of non-axisymmetric instabilities in rotating star will conclude these lecture notes. 



3.1 Linear Perturbations of Fluid Stars 

Fluid stars have two fundamental classes of oscillation modes that are referred to as nonra- 
dial and radial according to whether the perturbed quantities have an angular dependence or 
not. Furthermore, since not constrained by Birkhoff's theorem (cf. Section l2~Tj) . nonrotating 
(or spherical) stars are allowed to also have radial oscillations in addition to the nonradial 
ones. In the following sections I will discuss both radial and nonradial oscillations of nonro- 
tating stellar models, treating the first ones in full General Relativity, while restricting the 
discussion to Newtonian gravity for the second ones. 

As done for the perturbations of black holes, I first consider a fluid star which has a well 
defined "unperturbed" state which is one of equilibrium. The star will then be considered in 
a new "perturbed" state, not necessarily of equilibrium, which represents a small deviation 
from the initial one, i.e. 

f(t,r,e,<f>) = f (r) + f 1 (t,r,6, ( j>) , (108) 

where fo(r) is the value of the unperturbed quantity (considered here to be spherically 
symmetric) and f\ (t, r, 8, <fi) is its perturbation (considered here to be have also an angular 
dependence). Of course, linearity is guaranteed by the condition that /i//o <C 1- 

Two different approaches to the description of infinitesimal perturbations are possible 
and in which the physical variables are seen as 

• Eulerian perturbations: when the perturbed variables are considered at a fixed position 
in space with coordinates x, i.e. 



Sf(t,x) = f(t,x) - f (x) , 



(109) 
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where f(t,x) is the new value assumed by the physical variable (e.g. p, p) and fo(x) 
is the equilibrium value, independent of time. 

• Lagrangian perturbations: when the perturbation is obtained after comparing variables 
at two different spatial points occupied by the same fluid element 

Af(t,x) = f(t,x + i(t,x^)-fo(x ) , (110) 

where £(t,x) = x — xq is the Lagrangian displacement three- vector (£ = in the 
unperturbed configuration), f(t,x + £(t,x)) is the value of the perturbed physical 
variable, and fo(x ) the equilibrium value in the initial position x . 

A direct relation between the two descriptions is expressed through the following relation 
between the two operators A and S 

A = 5 + C (111) 

where £|* is the Lie derivative along £. 

3.2 Radial Oscillations: Relativistic Stars 

Consider the unperturbed star to be composed of a perfect fluid with four- velocity uq, energy 
density po and isotropic pressure p, whose stress-energy tensor takes the form 



ii 







Tp, = (po +Po)(u )^(u ) u +p Q g tlv , (112) 
where 

n = ^- , (113) 
m Q 

is the baryon number density of fluid particles of rest-mass mo and rest-mass density p*, so 
that e = Po/ P* — 1 is the unperturbed internal energy density. Note that, as in equation 
(J58|) . I have here used the upper zero index to indicate tensors referring to unperturbed 
quantities. The generic background spacetime of such a static spherical star is expressed 
through the line element 

ds* = - e 2uo ^dt 2 + e 2Ao(r W 2 + r 2 (d8 2 + sin 2 Odcj) 2 ) , (114) 

so that the Einstein equation in such a spacetime can be recast as a set of ordinary differential 
equations 

dpo Pornf 1 + ^\( 1 + ^\( 1 _2jny 1 ^ 



dr r 2 \ po J \ m J \ r J 



dv 1 dpo A , po 



-i 



, ,1 + - , (H6) 
dr po dr \ p J 

dm , 9 /,,„x 
— =4vrr 2 p , (117) 
dr 



27 



where the function m = m(r) has been introduced to re-express the g rr metric function, i.e. 

e 2A ° = (l - ^ 1 . (118) 

The function m(r) has the interpretation of gravitational "mass inside the radius r" because 
the integral form of f)117|) can be written as 

f R 

M = / A7ir 2 p dr , (119) 



where R is the radius of the star, so that M represents the total mass-energy of the star. 

The system of equations ()115j) - (jll7|) in the variables p , vo, p , and m is usually referred 
to as the set of Tolmann-Oppenheimer-Volkoff (or TOV) equations [2H1 121] • An additional 
equation is necessary to close the system and this is given by an equation relating the energy 
density and pressure appearing in the stress-energy tensor. Such a relation is usually referred 
to as equation of state (EOS) and in the case of a fluid in local thermodynamic equilibrium 
there always exists a relation of the type 

Po =Po(po,s ) , (120) 

where So is the specific entropy. If the change in the fluid entropy are very small so that 
the latter can be considered a small constant quantity, the EOS can be simply written as a 
relation between the energy density and the pressure and is referred to as barotropic. A very 
common example of a barotropic EOS is a polytropic EOS 

Po = Kpl , (121) 

where K and V are the polytropic constant and exponent, respectively. 

Once the "background" static stellar model has been constructed as a solution of the 

TOV equations, a radial perturbation of Eulerian type can be introduced in both the metric 
variables 

v Q (r) — > v(t,r) = u Q (r)+5u(t,r) , (122) 

A (r) — > A(t, r) = A (r) + 5X(t, r) , (123) 

and in the fluid variables 

Po(r) — > p(t,r) = p (r) + 5p(t,r) , (124) 

p (r) — > p(t,r)=p (r) + Sp(t,r) , (125) 

n o{ r ) — ► n(t, r) — no(r) + 8n(t, r) . (126) 



where I have assumed the perturbations to be isentropic, so that 5K = = ST. 
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Note that as a result of ()122JI - (|123|) . the metric line element ()114|) takes the non-diagonal 
form 

ds 2 = -e 2u{r '' l) dt 2 + lAdr'dt + e 2A(r '' t] dr' 2 + r' 2 (d6 2 + sin 2 6d(j) 2 ) , (127) 
which, however, can always be recast in a diagonal form of the type 

ds 2 = _ e 2u(t,r) dt 2 + e 2\{t,r) dr 2 + ^^2 + sin 2 ^2) ^ (128) 

after an appropriate coordinate transformation. 

Because only motions in the radial direction are possible, the problem has only one degree 
of freedom and it is convenient to use the radial component of the Lagrangian displacement 
vector = (£(£, r),0,0) to express all of the equations in terms of this scalar function. In 
particular, since the perturbed Eulerian velocity has components u r = + 5u r and u l = 
Uq + 5m*, the usual normalization conditions 

u^Ufi = -1 , = -1 , (129) 

together with the relation 

Au a = ^u a u b u c Ag bc , (130) 
can be used to derive the following expressions for the perturbed components 

5u r = e~ uo i , (131) 

Su f = e~ u °5v . (132) 

Expression (J131)) shows that, at first order in the perturbations, the time derivative of the 
Lagrangian radial displacement £ is related to the radial component of the three-velocity, 
i.e. 

u r Su r dr ■ , 
u l u dt 

where dr/dt is calculated along the worldline of the unperturbed star. The set of perturbative 
equations to be solved is therefore 

(5 [G^v — SirT^] = Einstein equations, (134) 

S [i^V a T™] = energy conservation, (135) 

S [P^V a T°^ = momentum conservation, (136) 

5 [V^nu 11 ) = 0] baryon number conservation, (137) 

— — = Ti — — adiabatic condition, (138) 

p n 
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where = (uo)^(u ) u + g^ v is the projection tensor orthogonal to u and 

dlnp 



ri = r x (r 



dlnp 



(139) 



and is assumed to be a constant in time. 

The full set of perturbative equations can be derived after expanding equations (|134j) - 
(I137j) and retaining only the first-order terms. In particular, the baryon conservation equa- 
tion ()137|) has a perturbed expression given by 



An = —n 



rV°0' + SX 



2 p Xq 



(140) 



with the prime indicating a partial radial derivative. Similarly, the perturbed energy and 
adiabatic conservation equations ()135|) and ()138|) yield the following relation between the 
perturbed energy density and baryon number 



a Po + Po A 
Ap = An . 

n 

Equation (j!41j) can also be expressed in terms of Eulerian perturbations to give 

1 



tip = -(Po+Po) 



8p = -Tipo 



rV°0' + SX 



2 p Xq 



rV°0' + 5X 



£Po- 



Next, the perturbations in the metric functions described by the field equations 
5X= -4n(p +p )re 2Xo Z , 

5u' = -^El e ^o + [Airp' r - An( Po + p )]e 2A ^ • 

Finally, the perturbed Euler equation ()136J) yields the perturbed "equation of motion' 



(141) 



(142) 

(143) 

yield 
(144) 

(145) 



(po + Po)e 2ix °- Uo) ii = -5p' - (5p + 5p)u' - (po + PoW , (146) 

whose right-hand-side shows the well-know restoring forces due to gradients in the pressure 
or in the gravitational potentials (both the unperturbed and the perturbed one). 
It is now convenient to introduce the following auxiliary variables 



rj = r e 



W = (Po+Po)- 



3Ao+^o 



> 



A = I>o- 



d Aq+3i^o 



> 



Q 



3 Ao+3^o 



(Pi 



o) 2 



4Pc 



D 2A 



.(po+Po)?" 2 r 



— - 8vr(p +p )po- 



(147) 
(148) 

(149) 

(150) 
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in terms of which equation (|146|) can be written [after using eqs. (jl42|) - (jl45j) ] as an inhomo- 
geneous wave equation of the type 

Wr} = (ATj'y + Qfj. (151) 

As done in Section |2~T1 for the perturbations of a black hole, a harmonic time dependence of 
the perturbations can now be assumed in terms of a constant complex frequency a 

rj(t,r) = r](r)e- ir7t , (152) 

so that equation (j!51|) becomes 

{An')' + (Q + a 2 W)n = . (153) 

Equation (|153|) represents a typical Sturm- Liouville eigenvalue equation and, together with 
suitable boundary conditions, its solution provides the eigenvalues and the eigenfunctions 
for the radial perturbations. Fortunately, for spherical stars the boundary conditions at the 
extrema of the range where the eigenvalue problem needs to be solved are rather straight- 
forward to formulate and can be summarized as follows: 

• when r — > 0, all quantities should be regular (i.e. either zero or finite) and this then 
implies that £ and £' should be regular or, equivalently, that n ~ r 3 for r — > 0. 

• when r — > R, on the other hand, the perturbation should leave the Lagrangian pressure 
perturbation equal to zero at the stellar surface, i.e. Ap = for r — > R. In practice, 
then, this implies that 

e v ° 

-—T lPo n' = for r -> R . (154) 

Despite the fact that determining the eigenfunctions and eigenfrequencies of radial oscil- 
lations represents the simplest possible eigenvalue problem, doing this in General Relativity 
requires the numerical solution of equation ()153|) . I will not discuss here the details of how 
this is done in practice, but simply summarize in what follows the main results [2*B1 IT| l2*Bj: 

• The eigenvalues a n form an infinite discrete sequence with integer index n, i.e. 

ctq < a\ < of < • • • • (155) 

• The eigenvalues a 2 are all real. 

• The eigenfunction corresponding to the lowest eigenfrequency n Q (r) has zero nodes for 
< r < R, while the n-th eigenfuction n n (r) has n nodes in the same interval. 

• The eigenfunctions n n (r) are orthonormal with weight function w(r), i.e. 

rR 

/ n n n m w(r)dr = 5 nm . (156) 
Jo 
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When using a description of the Sturm-Liouville problem in terms of a variational 
principle, the eigenvalues can then be calculated as 



a 2 = extremal of 



!*[A{j) 2 -Qn 2 ]dr 
wrfdr 



(157) 



While formally elegant, the variational method is not very used in numerical calcula- 
tions 



The absolute minimum of (|157|) signs the frequency of the fundamental mode of pul- 
sation. If this quantity is negative, the mode will grow in time on a timescale r g [i.e. 
T) ~ exp(— icrt) ~ exp(t/r g )] and the star will be unstable to radial oscillations. If, on 
the other hand, this quantity is positive, the mode will decay on a timescale Td [i.e. 
n ~ exp(— iat) ~ exp(— t/r d )] and the star will be stable. Since the denominator of 
(I157|) is positive-definite, the stability against radial perturbations is given by 



radial stability a 2 > / [A{n') 2 - Qn 2 ]dr > . (158) 

Jo 

It is also very instructive to consider the solution of equation (|157|) in the first-order 
post-Newtonian (1PN) approximation of General Relativity, in which all of the variables are 
expanded up to first order in the ratio M/R. In this case, the fundamental frequency of 
oscillation <t (not to be meant as an unperturbed quantity!) can be written as 

3\W\ - 

^ 2 = ^(ri-r c ), (159) 

where \W\ is the star's gravitational binding energy, I = An j (p Q r 2 )r 2 dr is the trace of the 
second moment of the mass distribution and fi is the pressure averaged adiabatic index 



Air J^Tip r 2 dr 



■>R- 

Tt EE - - . ( 160 ) 

An J Q p r 2 dr 

Equation ()159|) can also be interpreted as a criterion stating that the stability of the 

star to the fundamental mode of radial oscillation depends on the values of the averaged 
adiabatic index T\ when compared with the "critical" adiabatic index 

Tc = t + ' (161) 
where a > is a constant defined in terms of the equilibrium quantities 

Atx R f R ( ml m \ . 

a ee — — — / 3p ^ + 4p — r 2 dr , (162) 



3 M\W\ J \ r r 2 ^ r 

and mo(r) is the rest-mass inside the radius r. As a result, the stability of a relativistic 
stellar model to radial perturbation can be summarized as 

f i > T c => the star is stable to radial perturbations, 
fi < T c =^r- the star is unstable to radial perturbations, 
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Since for a Newtonian polytropic star the critical value of the adiabatic exponent discrim- 
inating stability is r c | Ncwt =4/3, equation (|159|) expresses the fact that the 1PN correction 
(and therefore the effects related to a stronger gravitational field) are 0(M/R), thus making 
the criteria of stability more severe (r c | Nowt < r c | 1PN ). Stellar models that are radiation- 
pressure dominated and very massive white dwarfs are very well described by polytropes 
with Ti = 4/3 and therefore the post-Newtonian contributions expressed in (j!59|) are very 
important to assess the stability of these objects against radial perturbations. 



3.3 Nonradial Oscillations: Newtonian Stars 

The discussion of the properties of nonradial oscillation modes is sufficiently complex (already 
for nonrotating stellar models) that a first approach is best done within a simpler scenario. 
For this reason, hereafter the discussion will be made in terms of a Newtonian description 
of physics. The interested reader will find a detailed treatment of this subject in [27] and a 
review on nonradial oscillations in relativistic stars in |2*5] . 

Consider therefore an equilibrium state of a nonrotating and zero temperature star that 
can be described by the distribution of pressure, density and gravitational potential, functions 
of the radial coordinate r only, i.e. 

Po = Po(r) , Po = Po(r) , % = %( r ) ■ (163) 

These functions are obtained as solutions of a system of differential equations describing the 
equilibrium 

1 dp d$ , 
— f ~ -j- = , (164) 
p dr dr 

V 2 $ = AnGp , (165) 



P = P{p) , (166) 

and in the case of a polytropic EOS, this system can be cast into a second order differential 
equation, the Lane-Emden equation, for p or p [29J. The solution of the system ()164j) - (jl66p is 
possible after suitable boundary conditions are imposed at the centre (where the regularity of 
all the quantities should be required) and at the surface of the star (where the gravitational 
potential should be matched to the vacuum solution behaving as r _1 ). Introducing now the 
Eulerian perturbations, the different variables can be expanded as 

p(t, r, 9, <P) = po(r) + 8p(t, r, 9, 0) , (167) 
p(t,r,9,<l>)=p (r) + 5p(t,r,e,<f>) , (168) 
v(t, r, 9, 0) = v (r) + 5v = 5v(t, r, 9, 0) , (169) 



$(t,r,9,(f)) = $ (r) + 5<S>(t,r,6,4>) , 



(170) 
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and it is useful to compare expressions (jl67|) - (jl70j) with the corresponding relativistic ex- 
pressions (|122|) - (jl26|) to appreciate the nonradial aspect of these perturbations. 

Using the Lagrangian displacement vector £, the Lagrangian velocity perturbation of the 
stellar matter is [30J 

Av l = d t C, (171) 
and in the nonrotating stellar case (vl = 0) it reduces to 

5v l = d t C. (172) 

Note that here the Lagrangian displacement does no longer need to have a radial component 
only [cf. equation ([184)1 ] . Introducing now the perturbed quantities 5p, 5p, 5$ and 5v l in 
equations (|164j) - (|166|) and in the equation for the conservation of baryon number, one obtains 

5p + V( P oO = 0, (173) 



-a 2 i+ —V(5p) + V(5$) + — V($ ) = , (174) 
Po Po 

^d r [r 2 d r (5$)] + Vi(5$) = AnG5p , (175) 

^ = . (176) 

V P 

where is the two-dimensional Laplacian operator in spherical polar coordinates (i.e. a 
two-dimensional gradient operator on the r = const. 2-sphere) 

Vi = — [sin 6d e (sm 6d e ) + <9j] . (177) 
r z sin 6 r 

Using now the relation between the Eulerian and Lagrangian perturbations (jlllj) . equation 
(|176|) can be rewritten as 

5p + £(dp /dr) _ T ^p + £(dp /dr) 
Po Po 

or, alternatively, as 

S S- = ^-M, (179) 
P r x p 

with 

_ cilnpo 1 d\np ( 1 1 \ d\np 

A = ^r-T l -^r = \f-Y l )^r ^ (180) 

and with the second equality in (jHQj) applying for a polytropic EOS only. 
The quantity A defined in equation (jl80p is the Schwarzschild discriminant and is used in 
discussing of the convective instability in a star [3T]. In particular, if A > somewhere 
in the star (or equivalently if Ti > T), the matter there will be unstable to convective 
motions. The Schwarzschild discriminant is therefore related to restoring buoyancy forces 
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and is important to determine the frequency of gravity waves in the radial direction or, 
equivalently, the frequency at which a fluid element may oscillate around its equilibrium 
position under the influence of gravity and pressure gradients. Such a frequency is called the 
Brunt- Vdisdld frequency and is defined as 

N 2 = -Ag , (181) 

where g is the local gravitational acceleration 

_ GM(r) I dp, 
g = — = — — . 182 

r z po dr 

A simplification should now be introduced to handle the non-trivial angular dependence 
in the perturbations ()167|) — f)17Uj) . To this scope I will assume that all of the relevant per- 
turbative equations can be rewritten after decomposing the variables into vector spherical 
harmonics. In particular, given a three-vector V with separable variables, it can be decom- 
posed into its "polar" and "axial" parts [cf. eq. (|58|i] as 

V(r, 6, 0) = ^[Vp!(r)y /m e r + V P2 (r) V ± Y lm ] + ^["W^r x V±Y lm )\ , (183) 

l,m l,m 

where the first two terms represent the "polar part" (with radial coefficients Vpi and V^), 
while the third one represents the "axial part" (with radial coefficient Va)- Here, e r is the 
radial unit vector and, as a result, the axial terms in ()183|) have only components in the eg 
and directions. 

A considerable simplification in the treatment of nonradial oscillations in Newtonian 
nonrotating fluid stars comes from the fact that no axial term appears in Euler equations to 
produce an axial restoring force. Stated differently, if we decomposed £ as 

i = £e r + £l = e- lat W)Y lm e r + ^(r)V ± Y lm ] , (184) 

then all of the terms of equations ()173|) - (jl75|) involving £j_ can be eliminated using the 
condition [cf. eq. ()174|) ] 

-a 2 U + V± f — + 5$ J = . (185) 



.Po 

As a result, only the polar parts need to be considered. (Note that this is not the case for 
stars with a solid crust, or with rotation, or with a magnetic field.). 

Assuming now the usual harmonic dependence in time, Sp, Sp, <5$,£ ~ e~ lat , and noting 
that equations with different quantum numbers l,m are decoupled (so that the sums can be 
replaced by a single pair of l,m indices) one obtains 

r A dr cj. \ a z J pcj, a z r z 

par pc A v ' dr 



r z dr \ dr J r z 



VPcf 9 ) 
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where c 2 



is the local sound speed 




(189) 



and Lf 



is the square of the Lamb frequency for the l-ih mode 




(190) 



or, equivalently, L~[ represents the typical crossing time of sound waves of wavelength ~ 
2nr jl on a r = const, surface. 

The set of equations ()186jl - ()188)l fully describes the nonradial oscillation modes of a non- 
rotating Newtonian star and its solution requires, in general, a numerical calculation. I will 
not discuss here the details of how to obtain such a numerical solution but will concentrate, 
in the following Section, on using a simplified form of equations (|186|) - (jl88|) to derive a first 
classification of the different oscillation modes that are solutions of this set of equations. 

3.4 Classification of Stellar Oscillation Modes 

The initial step in the description of a complex phenomenology always consists of a "classifi- 
cation" and stellar oscillation modes are no different in this respect. In the case of nonradial 
modes of nonrotating stellar models, the classification is particularly simple when the Eu- 
lerian perturbations of gravitational potential can be neglected (i.e. 5$ = 0), so that the 
perturbed Poisson equation ()188j) is trivially satisfied. This is usually referred to as the 
Cowling approximation |H2] and for a fluid configuration which is not self-gravitating (i.e. 
one whose gravitational potential is such that $o + 8$ — 0), the Cowling approximation 
is actually an exact description of the pulsations [33]. The error made in employing the 
Cowling approximation depends rather sensitively on the specific mode investigated and for 
modes with predominant axial nature, this error can be rather small (i.e. < 10%), while it 
can become more significant with predominantly polar modes [3*H I35j. 

As customary in the analysis of local oscillations, I now introduce two new variables x 
and ?7, replacing 5p and £, and defined as 




(191) 




(192) 



As a result, equations (|186|) and (|187|) can be rewritten in the matrix form 




V 



o 2 — N 2 
r 2 h 
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dr 



(193) 



where 

hf^~ \ r ( N2 9 

h(r) = exp / - 

Uo V 9 cf 

Next, I consider the perturbations to have a harmonic spatial dependence in the radial 
direction expressed as 

X,V~e** r , (194) 

where k r is the radial wavenumber, assumed to be larger than the length-scale of radial 
variations in the equilibrium configuration (this is also known as the WKB approximation 
and is a distinctive feature of a "local approach"). Introducing this ansatz in equation ()193|) . 
one obtains the following dispersion relation for nonradial oscillations in nonrotating stars 

k 2 = ^ 2 - Ll) (- 2 - N 2 ) • (195) 

If k 2 > 0, the wavenumber k r is real and two travelling waves will exist. The physical 
interpretation of the dispersion relation f)195j) is simpler in the two limiting cases in which 
either N = or L\ = 0, so that the dispersion relation assumes the simpler forms 

a 2 = L 2 + k 2 r c 2 s , (196) 

and 

a 2 = N 2 + k 2 r c 2 . (197) 

Furthermore, in the long wavelength limit k r = 0, the two frequencies in (|196|) - (|197|) corre- 
spond respectively to the Lamb frequency o = L\ [cf. eq. (|19U|) ] and to the Brunt- Vaisala 
frequency a = N [cf. eq. ()181j) ]. 

If, on the other hand, k 2 < 0, the wavenumber k r is imaginary, corresponding to waves 
that are exponentially damped somewhere in the star and are thus said to be "evanescent". 

All of this is summarized in Figure El which shows the frequencies of the different modes 
of oscillations as a function of the radial position within the star. Such a diagram is also 
referred to as the propagation diagram of the dispersion relation (j!95J) for a given number I. 
The thick solid lines represent the values of the to the Lamb and Brunt-Vaisala frequencies 
in the long wavelength limit k r = and divide the diagram in four different regions. In two 
of these regions the waves are evanescent and the corresponding frequencies are indicated 
with dashed lines. In the remaining two regions, (referred to as the upper and lower branch, 
respectively), the wavenumbers are real and the corresponding waves can be trapped to form 
an eigenmode (or global) oscillation. The frequencies of such modes are indicated with solid 
lines and the number of filled dots indicates the position of the nodes in the corresponding 
eigenfunctions. In what follows I will use Figure |B] to provide a brief description of the 
different modes that can be classified in an ordinary star, distinguishing those with pure 
polar properties from those with pure axial ones. 



3.4.1 Even-parity (Polar) Modes 

• p-modes - They fill the high-frequency branch of the propagation diagram and have 
pressure gradients as the main restoring force. The eigenfunctions of a mode with 
quantum number n has n nodes. As n increases, the frequency increases and the 
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r/R 1 



Figure 6: Schematic propagation diagram in a simple stellar model. The two thick lines represent the Lamb 
and Brunt- Vaisala, frequencies respectively, while the horizontal lines show the location of eigenfrequencies 
of several modes. Dashed lines are used to show the position within the star at which the waves become 
evanescent, while solid lines refer to eigenmode solutions; filled dots on these lines show the radial location 
of the nodes of the corresponding eigenfunctions. 



wavelength becomes smaller. In the limit of short wavelength, these modes represents 
simple acoustic waves travelling in the star. In this same limit the frequency of the 
mode will tend to infinity. 

For a "canonical" neutron star, i.e. a neutron star with mass M = 1.4M and radius 
R ~ 14 km, the typical frequency of the lowest order p-mode is a few kHz. 

• g-modes - They fill the low-frequency branch of the propagation diagram and have 
a buoyancy force (produced, for instance, by gradients in temperature, composition 
or density) as the main restoring force. The eigenfunctions of a mode with quantum 
number n has n nodes. As n increases, the frequency decreases and the wavelength 
becomes smaller. In the limit of short wavelength the frequency of the mode goes to 
zero. 

In neutron stars, two different types of g-modes are expected. One of them can be of 
modes trapped in the high density core region of the star and has a typical frequency 
of ~ 0.1 kHz. The other g-modes are expected to be trapped in the outer layers of the 
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star (i.e. in the fluid "ocean" above the crust) and have a typical frequency of ~ 10 
Hz. 

• f -modes - These modes have a character which is intermediate between those of p- 
and g-modes and are also referred to as the fundamental modes of oscillation. For 
each the frequency of this mode is between the lowest order g-mode (i.e. the highest 
frequency g- mode) and the lowest order p-mode (i.e. the lowest frequency p-mode). 
Note that for a given pair of quantum numbers (l,m), only one /-mode exists and its 
eigenfunctions have no nodes. For an incompressible homogeneous sphere, the /-mode 
represents the only oscillation mode and its frequency scales with the mean density 

1 21(1 + 1) ( M\ , , 

*« = V-srr U) ■ (198) 

The typical frequency of the lowest order /-mode for a canonical nonrotating neutron 
star is ~ 2 — 3 kHz with a decaying time of ~ 0.1 — 0.5 s. 



3.4.2 Odd-parity (Axial) Modes 

As mentioned in Section 13.31 this type of modes has zero frequency for nonrotating stars 
that do not have a crust or a magnetic field. When rotation is present, however, these modes 
are no longer degenerate and give rise to a rich and complex phenomenology. 

• inertial modes (including r-modes) - These modes have the Coriolis force as the main 
restoring force. In Newtonian stars, inertial modes have velocity fields described by 
a mixture of polar- and axial-parity components to first order in the slow-rotation 
expansion. The so-called "r-modes" are a special class of inertial modes and have 
axial-parity components only. The velocity eigenfunctions have very simple expressions 
and the perturbations in the density and pressure appear at orders higher than the 
first one in the slow-rotation approximation. The Newtonian r-modes are indeed a 
generalization of the Rossby waves, well-known in the geophysical hydrodynamics and 
follow a simple dispersion relation that, in a frame corotating with the star, takes the 
form 

°- = ~W+T) ■ (199) 

where Qi is the rotational angular frequency of the star as measured in an inertial 
frame. Since a r Q{ < 0, the wave pattern is retrograde in the rotating frame and in 
an inertial frame (i.e. in the reference frame of a distant observer) the mode will be 
observed to have frequency 

Oi = mVti - o~ r = ^~ 1 ^ Z + 2 ^ f . (200) 

A good starting point in the analysis of the classical r-modes in Newtonian stars can 
be found in [36 J, while a review of the impact of r-modes in the development of an 
instability can be found in jSTj. Finally, a discussion of the interaction of r-modes 
with the magnetic fields that are likely to be present in a neutron star can be found 
in 
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3.4.3 Purely Relativistic Modes 

As discussed in Section Q a general relativistic spacetime has its own dynamical degrees 
of freedom, thus allowing for the existence and propagation of gravitational waves. This 
distinctive feature of General Relativity gives rise to stellar oscillation modes that do not 
have a counterpart in Newtonian stars. These modes are called w-modes and were first shown 
to exist only rather recently |4U| . These modes are characterized by spacetime fluctuations 
that couple with the matter only very weakly. As a result, the stellar fluid is left almost 
unperturbed during w-mode oscillations. These purely relativistic modes can be further 
classified according to their properties which I briefly summarise below. More information 
can be found in the review work by Kokkotas and Schmidt [12] . 

• Curvature Modes: they are the standard w-modes |1U] and are reminiscent of the 
QNMs in black holes. The main features of these modes is that they have rather high 
frequencies of ~ 5 — 12 kHz and very rapidly decaying times of ~ 0.02 — 0.1 ms, with 
the damping rate increasing as the compactness of the star decreases. Modes of higher 
order have higher frequencies and shorter damping times. 

• Trapped Modes: they exist only for supercompact stars (R < 3M). Their existence 
is due to the fact that for such compact stars, the stellar surface is inside the peak 
of the gravitational field's potential barrier [IT] (cf. Fig. 0J). The trapped modes are 
finite in number (i.e. ~ 7) and their number increases as the potential well becomes 
deeper (i.e. as the star becomes more and more compact). Trapped modes differ from 
curvature modes in that their damping times are considerably longer (i.e. a few tenths 
of a second), with frequencies that are from a few hundred Hz to a few kHz. At present 
it is not yet clear whether such ultra-compact stars can be built through a realistic 
equation of state. 

• Interface Modes: they resemble acoustic waves scattered off a hard sphere but do not 
induce significant fluid motion [32] ■ For typical neutron stars they have frequencies in 
the range 2 — 15 kHz and are extremely rapidly damped, with damping times are of 
the order of less than a tenth of a millisecond. 

More on w-modes will be discussed in the following Section. 

3.4.4 Some Empirical Expressions for /-, p- and it; -modes 

As mentioned in the previous Sections, the different modes of oscillation depend on the 
physical properties of the stellar model and can provide detailed information on them. Nev- 
ertheless, the frequencies of the lowest order modes are not very sensitive on the fine details 
of the stellar structure and can be estimated through simple empirical formulas. This ap- 
proach was first suggested by Andersson and Kokkotas [Hi] and has lead to the different 
empirical expressions I will discuss below. 

Since the relation between the square of the /-mode frequencies and the mean density is 
linear [cf. equation (|198|) ]. introducing the following normalized quantities: M = Mj (1.4M ) 
and R = Rj (10 km), allows one to express the /-mode frequencies through the simple relation 




1/2 



kHz , 



(201) 
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which shows that the typical /-mode frequency is expected to be around 2.4 kHz. Equally 
important is to deduce a corresponding relation for the damping rate of the /-mode and this 
can be roughly estimated as 



(oscillation energy) 



(power emitted in gravitational waves) 



(202) 



where the power emitted in gravitational waves can be deduced from the quadrupole formula. 
Doing so leads to a relation of the type 



1 



22.85 



I 1.G5 ( 



,-i 



(203) 



which indicates that a typical /-mode oscillation is damped over about 0.6 s. 

A similar procedure can be followed for the oscillation frequency of £>-modes, which can 
be roughly estimated to scale as 



M 



1.75 + 5.59 



M 
R 



kHz 



(204) 



so that the typical p-mode frequency is around 7 kHz. Unfortunately, the damping of p- 
modes is sensitive to changes in the modal distribution inside the star and different EOSs 
lead to rather different p-mode damping rates. As a result, no robust empirical expression 
can be derived in this case. 

This is clearly not the case for w-modes that do not excite a significant fluid motion and 
are therefore only weakly dependent on the characteristics of the fluid. Rather, it can be 
shown analytically that the frequency of w-modes is inversely proportional to the size of the 
star 40J, while the damping time is a function of the compactness [44J. As a result, rather 
robust relations can be found for the frequency of the first w-mode 
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and for its damping time 



1 

M 



5.74 + 103 



M 
~R 



67.45 



M 
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Expressions ()205j) and ()206j) show that w-modes have typical frequencies around 12 kHz (but 
these can vary sensibly with the radius of the star) and a damping time of ~ 0.02 ms (which 
is comparable to that of an oscillating black hole with the same mass). 



3.5 Non-axisymmetric Instabilities 

Some of the non-axisymmetric modes of oscillation in rotating stars may not be damped, 
but have amplitudes that grow exponentially in time. When this is the case, the oscillations 
are said to be "unstable" and the resulting instability can either be dynamical, if it develops 
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on the timescale set by the rotation or by the free-fall, or secular, if it develops on a much 
longer timescale set, for instance, by dissipative processes. Dynamical instabilities differ 
considerably from secular ones in that they are purely hydrodynamical, while the latter 
are triggered by dissipative process such as viscous dissipation, emission of gravitational 
or electromagnetic radiation, thermal losses, etc.. In both cases, however, the instabilities 
reflect the attempt of the rotating star to find a lower energy state either by changing its 
mass distribution (e.g. through variations of the momentum of inertia) or by violating the 
conservation of some quantity (e.g. circulation or angular momentum). A quantity which is 
often used to measure how close the rotating star is to the onset of an instability is the so 
called rotational parameter (3 

_ (rotational kinetic energy) _ T 1 / Q 
(gravitational energy) \W\ 3 \fi K 

The last equality in ()207j) has been derived for a Newtonian star, with Q/Q K being the stellar 
angular velocity normalized to the Keplerian value, that is, the value of the angular velocity 
at which matter can be shed at the stellar equator. Indicating with p the average rest-mass 
density, the Keplerian angular velocity can be estimated to be Q K ~ (2/3)y / 7rp. 

The parametrization (|2()7jl is independent of the rotation law and is particularly use- 
ful for differentially rotating objects. By definition and invoking the Virial theorem, the 
parametrization is constrained to be between /3 — (for a spherical object) and (3=1/2 
(for an infinitely extended, thin disc at rest) [31] • A well-known application of the rotational 
parameter ()207|) is offered by classical result for the onset of the dynamical instability in 
Newtonian rotating stars, and which has been estimated to be T/|W| ~ 0.27 for a variety of 
different equations of state and rotation laws. 

Since the secular non-axisymmetric instabilities are triggered by dissipative mechanisms, 
their development will be different according to whether they are driven by viscous processes 
or by the emission of radiation (either gravitational or electromagnetic). When viscous dissi- 
pation processes are present and the radiative losses are negligible, an initially axisymmetric, 
incompressible rotating object, i.e. a Maclaurin spheroid, will be deformed into a Jacobi el- 
lipsoid, i.e. into a uniformly rotating, homogeneous configuration with ellipsoidal surfaces 
(they are similar to rotating american "footballs"). This happens at roughly T/|W| ~ 0.14 
and is referred to as the viscous-driven /-mode instability (see j3H] for a complete discussion). 

When viscous processes are negligible, on the other hand, the growth of the non-axisym- 
metric modes can be driven by the emission of gravitational or electromagnetic radiation 
(although the latter is usually much smaller than the former). The instability that develops 
in this way is the so called CFS (Chandrasekhar-Friedman-Schutz) instability |JBJ|3U] and is 
produced by the coupling between the loss of energy and angular momentum via radiation 
and the non-axisymmetric oscillations modified by the stellar rotation (A qualitative intro- 
duction to the CFS instability is presented in Appendix A). Also in this case, an initially 
axisymmetric Maclaurin spheroid will be deformed into a uniformly rotating, homogeneous 
configuration with ellipsoidal surfaces, called Dedekind ellipsoid. The difference between the 
Jacobi and Dedekind ellipsoids is that in the latter the ellipsoidal surfaces are supported by 
internal circulations but the shape is stationary as observed by an inertial observer (they 
are therefore similar to nonrotating american "footballs"). This happens again at roughly 
T/|W| ~ 0.14 and is referred to as the /-mode CFS-instability. 
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In practice, the development of non-axisymmetric instabilities is much more complicated 
than what discussed in the two limiting cases above, because both viscous and radiative 
losses are active at the same time in realistic stars. As a result, the modes that are driven 
unstable by viscosity and deform a Maclaurin spheroid into a Jacobi ellipsoid (i.e. the 
Jacobi modes) tend to be stabilized by the emission of gravitational waves (the star develops 
non-axisymmetric "ripples" to remove the excessive angular momentum via the emission of 
gravitational waves). At the same time, the modes that are driven unstable by the emission 
of gravitational waves and deform a Maclaurin spheroid into a Dedekind ellipsoid (i.e. the 
"Dedekind modes" ) tend to be stabilized by the viscous dissipative processes (the increased 
shear stresses, for example, tend to remove, with the aid of the shear viscosity, the non- 
axisymmetric "ripples" responsible for the emission of gravitational waves). Computing the 
delicate balance between these two mechanisms in regimes of strong gravitational fields, 
high-density matter and temperatures is extremely difficult and at the core of the present 
research on the emission of gravitational waves from instabilities in relativistic stars. 

Because of the very large amplitudes that the oscillation modes can reach when driven 
unstable, the amount of gravitational radiation emitted can become considerable and these 
unstable stars can then become promising sources of gravitational waves, potentially de- 
tectable by the gravitational- wave observatories now working or being under construction. 
A recent review on this fascinating area of research can be found in |47j . 
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A A Qualitative Introduction to the CFS Instability 

The instability was first discovered by Chandrasekhar j^O], and subsequently considered by 
Friedman and Schutz jHUj, who have shown its generic nature. While a formal proof of the 
criteria for the instability are rather involved [3U], qualitative arguments on the properties 
of the instability can be given simply using a couple of illustrative examples. 

Consider, therefore, a rotating star which is undergoing non-axisymmetric oscillations. 
For simplicity I will consider the simplest non-axisymmetric perturbation with mode numbers 
I = m = 2. Because the star is rotating, the properties of the perturbations can be considered 
both in a reference frame which is corotating with it (i.e. the "rotating" frame) or in a 
reference frame which is not rotating and is fixed with respect to, say, distant stars (i.e. the 
"inertial" frame). 



rotating frame 




Q r = G r /m<0 



Figure 7: Schematic view from the rotating frame. The unperturbed star is shown on the left and the 
non-axisymmetric perturbation on the right. f2 r = and oy/m are the angular velocities of the star and of 
the wave-pattern, respectively. Note that the star is effectively rotating as observed in an inertial frame and 
its structure, therefore, is deformed into a spheroid by the centrifugal stresses. 



I will first discuss what would be observed in the rotating frame; this is summarized 
in Fig. [7| where I have shown schematically the unperturbed star on the left and the non- 
axisymmetric I = m = 2 perturbation on the right. In the corotating frame the star has 
a zero angular velocity (i.e. Q r = 0), but is nevertheless deformed into a spheroid by the 
centrifugal force. The perturbation on the other hand, has a nonzero frequency o~ r and the 
corresponding m = 2 wave-pattern is seen to rotate with angular frequency oy/m which 
is, say, negative. By definition, such a mode is referred to as "retrograde" and because the 
perturbed star is rotating at an angular velocity smaller than the initial one, the mode has 
a negative angular momentum J < in the corotating frame. 
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The non-axisymmetric perturbation generates a time- variation of the stellar mass multi- 
poles (and/or of the mass-current multipoles) and the gravitational waves that are produced 
in this way, carry positive amounts energy at infinity. The angular momentum carried at 
infinity jgw , on the other hand, can either be positive or negative according to the sense 
in which the perturbation is seen to rotate in the inertial frame. I therefore need to con- 
sider how the perturbation is observed in the inertial frame which, I recall, is the frame in 
which quantities like the total amount of energy and angular momentum can be measured 
unambiguously. 

The "view" from the inertial frame is summarised in Fig. |H] where, again, I have shown 
schematically the unperturbed star on the left and the non-axisymmetric perturbation on 
the right. An observer in this frame will then see the star rotating at a nonzero, say positive, 
angular velocity fij > and the non-axisymmetric perturbation with a wave pattern that is 
also rotating with angular velocity Oijvn. 



Fi gure 8: Schematic view from the inertial frame. The unperturbed star is shown on the left and the 
non-axisymmetric perturbation on the right. Qi and Uijm are the angular velocities of the star and of the 
wave-pattern, respectively. Note that the star is seen to rotate and appears therefore as a spheroid. 

It is not difficult to realize that the direction in which the wave-pattern rotates depends 
on both f2j and ay through the relation between the frequencies in the two reference frames: 
(Tj = mQi — a r [cf. expression (|2UU|) ]. As a result, the wave-pattern can either be "dragged" 
forward (aj > > ay), or backward (a« < a r < 0) by the stellar rotation fV The obviously 
interesting case shown in Fig. |HJ is the one in which the mode is dragged forward (i.e. 
<Ji/m > 0) and the non-axisymmetric mode, which is then said to be "prograde", is seen to 
rotate in the same sense as the rotating star (Note that the "dragging" of the wave-pattern is a 
purely kinematical and Newtonian effect, fundamentally distinct from the general relativistic 
"dragging of reference frames".). When this happens, f2j<7 r < and the conditions for the 
onset of the CFS instability are met. In this case, in fact, and because the sign of the angular 




Q, >0 



C//m> 
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momentum lost is determined by the sense of rotation of the oscillation's wave-pattern, the 
prograde non-axisymmetric perturbation will carry to infinity positive amounts of angular 
momentum, i.e. j gw > 0. 

For the observer in the rotating frame, on the other hand, the total angular momentum 
of the mode J(t) = Jo — jgw(t) becomes increasingly negative because of the losses through 
jgw{t) that continuously reduce J(t), i.e. J(t) < Jq < 0. As a result, the initially small 
non-axisymmetric perturbation with negative angular momentum in the corotating frame, is 
driven to large amplitude oscillations with a progressively larger negative angular momentum. 
Such perturbation emits increasingly large amounts of gravitational waves, thus feeding the 
development of the instability. Using a pictorial analogue, the development of the CFS is 
similar to someone's debts that get larger as new expenses (with positive amounts of money) 
are made. The growth of the instability stops when either nonlinear or dissipative effects 
become important and transfer energy from the unstable mode into the other available 
channels. 
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